Entanglement Echoes in Quantum Computation 
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We study the stability of entanglement in a quantum computer implementing an efficient quantum 
algorithm, which simulates a quantum chaotic dynamics. For this purpose, we perform a forward- 
backward evolution of an initial state in which two qubits are in a maximally entangled Bell state. 
If the dynamics is reversed after an evolution time t r , there is an echo of the entanglement between 
these two qubits at time t e = 2t r . Perturbations attenuate the pairwise entanglement echo and 
generate entanglement between these two qubits and the other qubits of the quantum computer. 
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The development of new techniques which could en- 
hance the reliability of quantum computation is intrin- 
sically connected with the study of its stability. Ev- 
ery physical implementation of a quantum computer will 
have to deal with errors, due to the coupling with the 
environment or to an imperfect control of the computer 
hardware. Therefore an accurate study of the stability of 
a quantum computer, while it is running quantum algo- 
rithms, is demanded [lj. 

Entanglement is arguably the most peculiar feature of 
quantum systems, with no analog in classical mechan- 
ics. Furthermore, it is an important physical resource, 
which is at the basis of many quantum information pro- 
tocols, including quantum cryptography Q and telepor- 
tation Q . For any quantum algorithm operating on pure 
states, the presence of multipartite (many-qubit) entan- 
glement is necessary to achieve an exponential speedup 
over classical computation [jj. Therefore the ability to 
control entangled states is one of the basic requirements 
for constructing quantum computers. 

In this paper, we introduce a suitable method to char- 
acterize the stability of pairwise entanglement in quan- 
tum computation, by considering the echo of an initially 
maximally entangled pair of qubits. Namely, we as- 
sume that our quantum computer is initially in the state 
IV'o) = \®b) ® \x)- Here, the first two qubits are pre- 
pared in a maximally entangled Bell state while 
the other n q — 2 qubits are set in a pure state (|x)), and 
they are completely disentangled from the Bell pair (n q 
denotes the total number of qubits in the quantum com- 
puter). The state l^o) first evolves according to the given 
quantum algorithm, described by the unitary evolution 
operator IA. Then we invert the sequence of quantum 
gates that implement this algorithm, that is we apply 
W . In the ideal case we would reconstruct the initial 
state, since U^U\ipo) — IV'o)- However, due to noise and 
imperfections, the initial state \ipo) is not exactly recov- 
ered. In particular, the first two qubits are no longer in 
a Bell state, and therefore their pairwise entanglement 
echo is reduced. Conversely, this pair of qubits becomes 



entangled with the other qubits, thus generating multi- 
partite entanglement. In this paper, we study numeri- 
cally the attenuation of the pairwise entanglement echo 
in a quantum computer implementing an efficient quan- 
tum algorithm which simulates quantum chaotic dynam- 
ics. We point out that the entanglement echo simulations 
discussed in the following are close in spirit to the spin 
echo experiments in many-body quantum systems in the 
presence of perturbations Q . 

We study the entanglement echo for the quantum al- 
gorithm simulating the sawtooth map dynamics . The 
sawtooth map is a periodically driven dynamical system, 
described by the Hamiltonian 

H(9,n,r) = ^-^l £ S(r - jT), (1) 

j=-oo 

where (n, 0) are conjugated action-angle variables (0 < 
< 2ir). The time evolution r — » r + T of this system is 
classically described by the map n = n + k{9 — 7t), and 
6 = 6 + Tn, where the bars denote the variables after one 
map iteration. By rescaling n — ► p = Tn, one can see 
that classical dynamics depends only on the parameter 
K = kT. The classical motion is stable for — 4 < K < 
and completely chaotic for K < —4 and K > 0. The 
quantum evolution in one map iteration is described by 
the unitary operator U : 

|$ = U\yj) = e -* T " 2 /2 e lfe ^ 2 / 2 |V>, (2) 

where n = —id/dd (we set h = 1). The classical limit is 
obtained by taking k — > oo and T — > 0, keeping K = kT 
constant. We study map J2J on the torus < < 2tt, 
— 7r < p < iv. With a n q qubits quantum computer, we 
can simulate the quantum dynamics of the sawtooth map 
with N = 2 n " levels, and we set T = 2vr/iV. The effective 
Planck's constant of the quantum system is H e g ~ 1/N 
and the classical limit corresponds to n q — > oo (heff — > 0) 
2\. We focus on the case K = 5, which corresponds to 
the chaotic regime. 
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It is convenient to simulate map J2Jl by means of the 
forward-backward Fourier transform between 9 and n 
representations. While the classical fast Fourier trans- 
form requires 0{N log 2 N) operations, an efficient quan- 
tum algorithm has been found 0, which uses the quan- 
tum Fourier transform and simulates J3J in 0(n} q — 
(log 2 N) 2 ) elementary quantum gates per map iteration. 
Moreover, all n q qubits are used in an optimal way, that 
is no extra work space qubits are required. In this way in- 
teresting physical phenomena, like dynamical localization 
, cantori localization, and anomalous diffusion could be 
simulated already with less than 10 qubits. Therefore the 
quantum sawtooth map represents an interesting testing 
ground for quantum computation, and it is important 
to understand the limits to the quantum computation of 
this model due to noise and imperfections. 

To compute the entanglement echo, we start from the 
initial state 




|V>o> = |*u> ® |x> = -j= (|oo) 



|11))<8>|00...0> (3) 



and we perform a forward evolution of the quantum saw- 
tooth map J2J up to time t = t r , that is U — U tr (the 
discrete time t = t/T denotes the number of map itera- 
tions) . Then we compute the time reversal evolution up 
to the echo time t e = 2t r (namely, W — (£A) tr ). Our al- 
gorithm can be decomposed into single-qubit Hadamard 
gates and two-qubit controlled-phase shift gates Q{. In 
particular, the Hadamard gate can be written as n H ■ cr, 
where n H = (l/v2> 0, l/v2), and cr = (a x ,a y ,a z ), the 
CTj's being the Pauli matrices. Due to the imperfect con- 
trol of the quantum system during the quantum compu- 
tation, the initial state is not perfectly recovered. In this 
paper, we only deal with unitary errors, modeled by noisy 
gates. We assume that errors tilt the rotation axis h H 
by an angle randomly fluctuating in the interval [— e, e]. 
In the noisy controlled-phase shift gates, random phases 
of amplitude inside the interval [— e,e] are added. We as- 
sume that the errors affecting two consecutive quantum 
gates are completely uncorrelated. 

Since we consider unitary errors, the state \ip(t)) of 
the quantum computer at any time t is still a pure state. 
By tracing \ip(t)) over all the qubits, except those ini- 
tially prepared in a Bell state, we obtain the reduced 
density matrix p 12 (t) = Tr 3 ,..., n<! (\ip(t))(ip(t)\). Note 
that, in general, qubits 1 and 2 are no longer disentan- 
gled from the other qubits of the quantum computer, 
and therefore p± 2 is a mixed state. We evaluate the en- 
tanglement of formation E(t) of the state p±2 following 
Ref. [8| . First of all we compute the concurrence, defined 
as C — max(Ai — A2 — A3 — A4, 0), where the Aj's are the 
square roots of the eigenvalues of the matrix R = P12P12, 
in decreasing order. Here p\2 is the spin flipped matrix 
of pi2, and it is defined by pi2 = {o~ y ® o~ y ) p* 2 io y ® a y ) 
(note that the complex conjugate is taken in the com- 
putational basis {|00), |01), |10), |11)}). Once the con- 



FIG. 1: Entanglement echo in a noisy quantum computer im- 
plementing the sawtooth map algorithm in the chaotic regime 
at K — 5, with n q — 5 qubits (dashed line) and n q = 8 qubits 
(solid line), and perturbation strength e = 1CP 2 . We start 
from the initial state , and, from t = to t r — 20 a forward 
evolution of the sawtooth map is applied. After that time, we 
invert the dynamics. The echo occurs at time t e = 2t r — 40. 
Top: entanglement of qubits 1 and 2. Bottom: entanglement 
between these two qubits and the other ones. 



currence has been computed, entanglement is obtained 
as E = h((l + vl — C 2 )/2), where h is the binary en- 
tropy function: h{x) = — xlog 2 a; — (1 — a;)log 2 (l — x). 
We also compute the Von Neumann entropy S(t) = 
— Tr [/012(f) log 2 /012(f)] of the reduced density matrix P12. 
This quantity measures the entanglement between the 
qubits 1 and 2 and the other n q — 2 qubits of the quantum 
computer. In particular, we compute the entanglement 
echo E(t e ) and the Von Neumann entropy S(t e ) at the 
echo time t e . 

A typical numerical simulation of the entanglement 
echo is shown in Fig. ^ The upper part shows the be- 
havior of the pairwise entanglement of the two qubits 
initially prepared in a Bell state; in the lower part we 
plot the Von Neumann entropy of this two-qubit sub- 
system. The dynamics completely destroys the initial 
pairwise entanglement, which is partially recovered only 
at the echo time t e |9|- Instead, the Von Neumann en- 
tropy quickly reaches the saturation value. These re- 
sults can be understood as follows: in few map itera- 
tions (f ~ 1), the chaotic dynamics transforms the ini- 
tial state 1-00) hito an ergodic state. We can expand 
the state \tp{t)) at time t over the computational ba- 
sis: \lp(t)) = J2 ai ...a 

(f)|ai . . . a n ), where 
a, =0,1 for % = 1, ...,n q . Since \i/j(t)) is ergodic, the 
coefficients c a . 



(t) have random phases and ampli- 
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tudes \ca u ...,an I ~ 1/VN (to assure wave function nor- 
malization). For an ergodic system, the reduced den- 
sity matrix of a two-qubit subsystem is essentially diag- 
onal. Indeed, the diagonal matrix elements are given by 

(pi2)ai,a 3 ;oi,Qi2 = J2a 3 a „ \ c ai,...,a nq | 2 , and their value 

is S3 1/4, since they are given by the sum of N/4 positive 
terms, whose value is ~ 1/JV. The off-diagonal matrix 
elements of /912(f) are instead given by the sum of N/4 
terms of amplitude 1/N and random phases. Hence their 
value is 0(l/y/N). For such a nearly diagonal density 
matrix /912(f), the entanglement of formation can be an- 
alytically computed and we find that E — 0. This means 
that chaotic dynamics quickly destroys the entanglement 
of any two-qubit subsystem, as shown in Fig. Under 
the hypothesis that the wave function is ergodic, it is 
also possible to compute analytically the Von Neumann 
entropy S of the qubits 1 and 2. After averaging over 
noise realizations, one has S S3 2 — 8/(iVln2) [lO|,lll|, in 
good agreement with our numerical data. This value is 
close to the maximum possible entropy of the two-qubit 
subsystem, S m ax = 2. It is interesting to note that it is 
possible to invert the quantum dynamics after long times 
(much longer than the times for relaxation to statistical 
equilibrium) and recover the initial out of equilibrium 
state. This is a clear demonstration of the stability of 
the quantum motion in contrast to the high instability of 
the classical chaotic motion )l2j. 

We now focus on the stability of the entanglement echo 
at time f e . The numerically computed entanglement echo 
and Von Neumann entropy at the echo time f e are shown, 
for different noise strengths, in Figs.|2and|21 respectively. 
Fig-Elshows that noisy gates attenuate the entanglement 
echo, and, if the time t e is long enough, completely de- 
stroy it. This can be explained by noticing that unitary 
errors transform the echo state into a state which be- 
comes closer to an ergodic state, as t e increases. As we 
have seen above, an ergodic state of the whole quantum 
computer implies a pairwise entanglement E — [l3j |. 

The decay of the entanglement echo can be understood 
by considering that each noisy gate transfers a probabil- 
ity of order e 2 from the ideal state to all other states. 
Since there are no correlations between consecutive noisy 
gates, the population of the initial state decays exponen- 
tially, and we can write the echo state as follows: 

|V(te)) « e- C£2 "»^Vo) + £ a a (t e )\ a ), (4) 

where C is a constant to be determined numerically, n g t e 
is the total number of gates required to perform the echo 
experiment (n g = 3n 2 + n q being the number of gates 
per map iteration), and the sum runs over all the states 
\a) = \<xx...a n ) of the computational basis, except for 
the two states involved in the initial wave vector (^o) 
(\a A ) = |000...0) and \a B ) = |110...0)). Given the com- 
plexity of the dynamics simulated by the quantum algo- 
rithm for the sawtooth map, it is reasonable to assume 
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FIG. 2: Attenuation of the entanglement echo of a Bell pair 
in the sawtooth map, at K — 5, n q = 7, and, from right to 
left, e = 7.5 x 10" 3 , 10~ 2 , 1.2 x 10~ 2 , 1.5 x 1(T 2 , 2 x 10~ 2 , 3 x 
10 _2 ,4 x 10 -2 . Here and in the following figures data are 
averaged over 400 runs with different noise realizations. Inset: 
semilogarithmic plot of the same curves. 
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FIG. 3: Von Neumann entropy of the two initially entangled 
qubits as a function of the echo time f e , with same parame- 
ter values as in the previous figure. Inset: approach to the 
saturation value Soo for the same curves. 



that the coefficients a a have random signs and ampli- 
tudes of the order of - e- ^ 11 ^) / (N - 2) (to as- 
sure that |^(f e )) has unit norm). We can compute the 
entanglement echo E(t e ) from the expression (@J for the 
echo wave function. For e 2 n g t e -C 1, it turns out that 
E(t e ) w 1 - (3/21n2)Ce 2 n g f e . Therefore the entangle- 
ment echo is stable up to time f e oc l/(e 2 n g ) oc l/(e 2 n 2 ). 
This theoretical estimate is confirmed by Fig. in which 
we plot the characteristic time scale t* for the decay of the 
entanglement echo, defined by the condition E(t*) = c 
(we take c = 0.9). It is clearly seen that f* oc n~ 2 e~ 2 . 
Therefore noisy gates degrade the entanglement echo af- 
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of the quantum computer. We point out that, since the 
entanglement can be measured experimentally in an effi- 
cient way ^3 , entanglement echo experiments analogous 
to the numerical simulations discussed in this paper could 
be implemented in quantum processors with a small num- 
ber of qubits (4-10) and a few hundreds of gates. These 
experiments are close to present capabilities 0, 0, 0] 
and would bring new insights in our understanding of the 
limits to quantum computation due to decoherence and 
imperfections. 
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bility in quantum information processing" . 



FIG. 4: Time scale i* for the decay of the entanglement 
echo in the sawtooth map at K = 5, for different strengths 
and number of qubits: n q — 4 (empty circles), 5 (filled cir- 
cles), 6 (empty squares), 7 (filled squares), 8 (empty trian- 
gles), 9 (filled triangles), and 10 (diamonds). Straight line: 
t* — A/n^e 2 , with the fitting constant A m 6.04 x 10 -2 . Inset: 
rate V of the approach to equilibrium of the Von Neumann 
two-qubit reduced entropy. The straight line gives T = _Be 2 n 2 , 
with B w 2.34. Logarithms are decimal. 

ter a number n* e of elementary gates which is independent 
of the number of qubits (n* — n g t* cx e 2 ). 

The Von Neumann entropy S(t e ) of the reduced two- 
qubit subsystem at the echo time is shown in Fig. [3J It 
saturates, for sufficiently long echo times, to the value 
Soo « 2 — 8/(N In 2), as expected for an ergodic state of 
the quantum computer. We can compute the approach 
to equilibrium from Eq. that gives Soo — >S(i e ) oc 
exp(— 2Ce 2 n g t). This theoretical prediction is borne out 
by our numerical data shown in Figs. 13141 The inset of 
Fig.[3]shows that the approach to the asymptotic value is 
exponential. Indeed, we have S(t e ) ~ <S' 00 (1 — e~ rtc ). In 
Fig. 21 (inset) we plot the rate T for different number of 
qubits and perturbation strengths. We see that T cx e 2 n 2 . 

It is interesting to compare the entanglement echo de- 
cay with the decay of the fidelity, which is the usual 
tool used to characterize the stability of quantum com- 
putation P, Q ■ The fidelity / at time t e is defined as 
f(t e ) = \(ip(t e Mo)\ 2 , and Eq. implies that f(t e ) « 
exp(— Ce 2 n g t e ). This is in agreement with our numeri- 
cal data (not shown here). Therefore the decay of the 
fidelity, the entanglement echo and the approach to equi- 
librium for the reduced Von Neumann entropy take place 
in the same time scale oc l/(e 2 n g ) [lij . 

In summary, we have proposed a suitable method, the 
entanglement echo, to study the stability of entanglement 
under perturbations. We have shown that noise destroys 
the entanglement of a pair of qubits and produces entan- 
glement between these two qubits and the other qubits 
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